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ABSTRACT 

QSO absorption spectra provide an extremely useful probe of possible cosmological 
variation in various physical constants. Comparison of Hi 21 cm absorption with cor- 
responding molecular (rotational) absorption spectra allows us to constrain variation 
in y = a 2 g p where a is the fine structure constant and g p is the proton g-factor. We 
analyse spectra of two QSOs, PKS 1413+135 and TXS 0218+357, and derive values 
of Ay/y at absorption redshifts of z — 0.2467 and 0.6847 by simultaneous fitting 
of the Hi 21 cm and molecular lines. We find Ay/y = (—0.20 ± 0.44) x 10 -5 and 
Ay/y = (—0.16 ± 0.54) x 10~ 5 respectively, indicating an insignificantly smaller y in 
the past. We compare our results with other recent constraints from the same two 
QSOs (Drinkwater et al. 1998; Carilli et al. 2000) and with our recent optical con- 
straints which indicated a smaller a at higher redshifts. 
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1 INTRODUCTION 

Currently, many modern unified theories provide strong mo- 
tivation for experimental searches for variation in physical 
constants. Kaluza-Klein, Superstring and M-theory require 
the existence of extra spatial dimensions which are com- 
pactified on small scales. The coupling constants in our 3- 
dimensional subspace of these theories are related to the 
scale sizes of the extra dimensions. Therefore, if the ex- 
tra dimensions evolve with cosmology, we expect our 3- 
dimensional constants to vary (e.g. Forgacs & Horvath 1979; 
Marciano 1984; Barrow 1987; Damour & Polyakov 1994). 
Currently, there seems to be no mechanism for keeping the 
sizes of the extra dimensions constant (Li & Gott 1998) and 
so these modern theories naturally predict variations in the 
fundamental constants. 

Experimentally, quasar (QSO) absorption lines provide 
an ideal probe of cosmological variation of physical con- 
stants. Savedoff (1956) first analysed doublet separations 
seen in galaxy emission spectra to obtain constraints on vari- 



ation in the fine structure constant, a = e 2 /he. Absorption 
lines in intervening clouds along the line of sight to QSOs 
are substantially narrower than intrinsic emission lines and 
so provide a much more precise probe of a at high red- 
shift. Bahcall, Sargent & Schmidt (1967) first used the dou- 
blet spacings of gas seen in absorption against background 
QSOs to derive firm upper limits on possible variation of a. 
We summarize results from optical absorption line studies in 
Section suffice it to say here that the tightest constraint 
on Act /a found to date is a tentative detection of variation 
at Aq/q = (-0.72 ± 0.18) x 10~ 5 in 49 absorption clouds 
over the redshift range 0.5 < z < 3.5 (Murphy et al. 2001a; 
Webb et al. 2001). [] 

Even tighter constraints can be placed on variability 
of combinations of physical constants using radio QSO ab- 
sorption spectra. A comparison of Hi 21 cm and molecular 
rotational absorption spectra is particularly convenient. The 
Hi 21 cm hyperfine transition frequency is proportional to 
HpUb / {"ho?) where fi p = g v e%/ '(4m p c) and /is = eh/(2m e c) 



E-mail: mim@phys.unsw.edu.au (MTM) 



t See Murphy et al. (2001a) for a summary of all relevant con- 
straints on Aa/a. 
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is the Bohr magneton. Here, g p is the proton g-factor, m p 
and m e are the masses of the proton and electron respec- 
tively and a = h 2 /(m e e 2 ) is the Bohr radius. The rota- 
tional transition frequencies of diatomic molecules, such as 
CO, are proportional to h/(Ma 2 ) where M is the reduced 
mass. Therefore, the ratio of the hyperfine and molecular 
rotational frequencies is proportional to a 2 g p M/m p . Vari- 
ations in this quantity will be dominated by variations in 
y = a 2 g p since variations in M/m p are suppressed by a fac- 
tor m p /U ~ 100 where U is the binding energy of nucleons 
in nuclei. If any variation in y occurs, it will be observed 
as a difference in redshift between the Hi 21 cm (zh), and 
molecular (z mo i) absorption lines: 

Ay _ y z - y ^ Az _ z mol - z H 
y yo 1 + z 1 + Zmol 

where y z and yo are the values of y at the absorption redshift 
z and in the laboratory respectively. 

The first comparison of H 1 21 cm and molecular absorp- 
tion was made by Varshalovich & Potekhin (1996). They 
compared the published redshifts of the CO absorption (re- 
ported in Wiklind & Combes 1994) and Hi 21 cm absorp- 
tion (reported in Carilli, Perlman & Stocke 1992) towards 
PKS 1413+135. They interpreted any shift as a change in the 
molecular mass, M. Wiklind & Combes (1997) carried out 
a similar analysis. Drinkwater et al. (1998, hereafter D98) 
pointed out that such a shift actually constrains Ay/y rather 
than AM/M with a similar argument to that above. Var- 
shalovich & Potekhin's (1996) value of AM/M translates to 
Ay/y= (-4 ±6) x 10" 5 . 

Like Varshalovich & Potekhin (1996), in D98 we used 
the published redshifts for the CO(0— 1) lines towards 

PKS 1413+1350. However, we improved the precision of 
Ay/y by fitting Voigt profiles to the Hi 21 cm data rather 
than using the published redshifts. The errors on the redshift 
of each velocity component in the H I data were reduced by 
an order of magnitude compared to the estimates of Var- 
shalovich & Potekhin (1996). We also applied this analysis 
to another absorber towards TXS 0218+357 at z a b s = 0.685. 
In D98 we derived an upper limit on any variation in y for 
both absorbers: \Ay/y\ < 0.5 x 10" 5 . 

The present work aims at turning our upper limit on 
Ay/y in D98 into a measurement of Ay/y by fitting simulta- 
neously the H 1 21 cm data and the rotational lines of several 
molecular species. We describe the available data in Section 
2 and describe our analysis and results in Section 3. We 
discuss our new results in Section 4, comparing them with 
other constraints on a variability. 



2 AVAILABLE DATA 

Only 4 QSOs have had mm-band rotational molecular ab- 
sorption detected along their lines of site: TXS 0218+357, 
PKS 1413+135, B3 1504+377 and PKS 1830-211. Below we 
describe the molecular and H I data available for the first 

t The CO data of Wiklind & Combes (1994), used by Var- 
shalovich & Potekhin (1996), showed a — llkms -1 offset from 
the H 1 21 cm velocity. We used a new, corrected CO data set re- 
ported in Wiklind & Combes (1997) which shows no such offset. 



two of these QSOs. We do not have sufficient data to con- 
sider PKS 1830— 211. The molecular absorption profile is 
very broad (FWHM ~ 40kms"\ Wiklind & Combes 1998) 
and so constraints on Ay/y would be weak. We do not con- 
sider B3 1504+377 any further due to a 15kms _1 offset be- 
tween the Hi 21 cm (Carilli et al. 1997, 1998) and HCO+ 
(Wiklind & Combes 1996) velocity scales as noted in D98. 
Checks on the origin of this difference have not revealed any 
instrumental or human error and so it may well be that 
the molecular and H I absorption occurs in different clouds. 
However, as noted in D98, this seems unlikely because the 
molecular and Hi velocity structures are very similar. Fur- 
ther data is required to clarify this problem. 



2.1 TXS 0218+357 

TXS 0218+357 is a gravitationally lensed QSO (Patnaik et 
al. 1993), probably at a redshift z cm « 0.94 (Wiklind & 
Combes 1995), showing absorption in the lensing galaxy at 
Zabs = 0.6847 (Carilli, Rupen & Yanny 1993). VLBI ob- 
servations (Patnaik et al. 1993; Biggs, Browne & Wilkin- 
son 2001) show two, compact, flat-spectrum components (A 
to the SW and B to the NE) and a steep-spectrum Ein- 
stein ring. The background QSO shows intensity variabil- 
ity on intraday (Biggs et al. 2001) and longer (~monthly) 
time-scales (O'Dea et al. 1992; Patnaik et al. 1993). We use 
the spectrum of this QSO despite the fact that it is lensed 
since the A component provides the dominant absorption in 
both the molecular (Menten & Reid 1996) and 21 cm (Car- 
illi et al. 2000, hereafter COO) bands. Thus, the profile of the 
z = 0.6847 absorber should not vary significantly over the 
time between molecular and Hi 21 cm observations. 

In this paper we use the H 1 21 cm spectrum of the z — 
0.6847 absorption system published by Carilli et al. (1993) 
(FWHM = 6.9 kms" 1 ). Spectra of CO(l-2), 13 CO(l-2), 
C ls O(l-2) and CO(2-3) are taken from Combes & Wiklind 
(1995) (FWHM ~ 30ms" 1 ) and those of HCO + (l-2) and 
HCN(l-2) are from Wiklind & Combes (1995) (FWHM ~ 
2kms~ ). These data are presented in Fig. 1. 



2.2 PKS 1413+135 

PKS 1413+135 is at the centre of an edge on spiral galaxy 
at z om = 0.24671 and the absorption occurs in the disk of 
this galaxy (Wiklind & Combes 1997). VLBI observations 
(COO) show that the 21 cm continuum flux is dominated by 
a jet extending from an inverted-spectrum nucleus. The nu- 
cleus dominates the mm continuum emission (Perlman et 
al. 1996). Thus, the Hi 21 cm and molecular absorption lie 
along different sight lines separated by ~ Of! 03 (COO). We 
address this problem and its effect on our results in Section 

Here we use the Hi 21 cm spectrum of the z — 
0.2467 absorption system observed by Carilli et al. (1992) 
(FWHM PS 1.3 kms -1 ). The molecular absorption lines 
of CO(0-1), HCO + (l-2) and HCO+(2-3) (FWHM ~ 
40 ms -1 ) were observed with the narrow band autocorrela- 
tor at the IRAM 30 m telescope on Pico Veleta. The details 
of these observations are reported in Wiklind & Combes 
(1997). These data are presented in Fig. 2. 
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Figure 1. Spectra of the z = 0.6847 absorption system towards TXS 0218+357. The data are plotted as a histogram and the solid line 
represents our Voigt profile fit. The tick-marks above the data show the position of the fitted velocity components. The residuals (i.e. 
[data] -[fit]), normalized to a constant lc-error (defined by the signal-to-noise ratio in the continuum), are shown above each spectrum. 
Note that two independent spectra of CO(l— 2) have been included. 



3 ANALYSIS, RESULTS AND ERRORS 

We used vpfit^] to fit multi- velocity component Voigt pro- 
files to both the Hi 21 cm and molecular absorption lines. 
The laboratory values for the transition frequencies were 
taken from Pickett et al. (1998)0 and Essen et al. (1971) and 
are given in Table |l[ By simultaneously fitting several molec- 
ular transitions of differing strengths in the same absorber, 
we get the best possible estimate of the velocity structure. 
We determine the number of velocity components to be fit- 
ted to each system by requiring that the reduced x 2 for the 
fit be ~ 1. 

We use an iterative technique to find the best fit value 
of Ay/y. We vary Ay/y outside vpfit (i.e. we vary the ratio 
of the hyperfine and molecular frequencies) and fit all lines 
simultaneously. During the fit we force the redshift param- 



Sce 



http://www.ast.cam.ac.uk/~rfc/vpfit.html for details 



about obtaining VPFIT. 

If See http://spec.jpl.nasa.gov 



eters for corresponding velocity components to be the same 
for all transitions, thereby reducing the number of free pa- 
rameters. We change Ay/y at each iteration of the routine, 
i, to find Xi as a function of Ay/y. The best fit value of 
Ay/y is that which gives the minimum \ 2 \ Xmim an d the la 
error in Ay/y is found using Xmin + !• 



For the z = 0.6847 absorption system towards 
TXS 0218+357 we find the best fit value of Ay/y = (-0.16+ 
0.36) x 10 -5 and for the z = 0.2467 absorption system to- 
wards PKS 1413+135 we find Ay/y = (-0.20 + 0.20) x 10~ 5 . 
To check internal consistency, we have investigated the ef- 
fect of removing individual transitions from the analysis of 
each absorption system. Our values of Ay/y showed no sig- 
nificant change upon removal of any molecular transition. 
In particular, we found no change when we removed either 
or both the CS(2— 3) and 13 CO(0— 1) transitions from our 
analysis of TXS 0218+357. The low velocity component does 
not seem to be present in these transitions and was not fit- 
ted (see Fig. 1). The fact that the results are insensitive to 
removal of these lines indicates the robustness of our line 
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Figure 2. Spectra of the z = 0.2467 absorption system towards PKS 1413+135. The data are plotted as a histogram and the solid line 
represents our Voigt profile fit. The tick-marks above the data show the position of the fitted velocity components. The residuals (i.e. 
[data] -[fit]), normalized to a constant lcr-error (defined by the signal-to-noise ratio in the continuum), arc shown above each spectrum. 



Table 1. Laboratory frequencies of transitions used in our anal- 
ysis. Molecular frequencies are taken from Pickett et al. (1998) 
and the Hi 21 cm frequency is from Essen et al. (1971). 



Molecule 


Transition 


Frequency /GHz 


CO 


0-1 


115.2712018(5) 




1-2 


230.5380000(5) 




2-3 


345.7959899(5) 


13 co 


1-2 


220.3986765(53) 


c 18 o 


1-2 


219.5603568(81) 


cs 


2-3 


146.969033(50) 


HCO+ 


1-2 


178.375065(50) 




2-3 


267.557619(10) 


HCN 


1-2 


177.261110(2) 


Hi 


21cm 


1.420405751766(1) 



different sight lines with the molecular and Hi 21 cm ob- 
servations may result in an additional random error. D98 
have investigated this problem empirically by fitting several 
QSO spectra with Galactic absorption and comparing the 
fitted redshifts of the HCO + and Hi velocity components. 
D98 find a close correspondence between the velocity com- 
ponents of HCO + and H I with a Gaussian dispersion of only 
1.2 kms -1 . This corresponds to an error Ay/y = 0.4 x 10 -5 . 
Our two results, being of similar value and so close to zero, 
suggest that this error estimate may be too large. However, 
it is clear that such a conclusion can only be reached with a 
larger sample of absorption systems and so we have added 
this error in quadrature to obtain our final results: 



&V/V = (-0.20 ± 0.44) x 10" 5 at z = 0.2467 and 
Ay/y = (-0.16 ± 0.54) x 10 -5 at z = 0.6847 



(2) 



fitting method and highlights the advantage of using many 
different molecular transitions. 

The la errors we quote above are statistical only. How- 
ever, as noted in Section 2.2, the fact that we probe slightly 



Aft/a = (-0.10 ± 0.22) x 10~ 5 
Aa/a = (-0.08 ± 0.27) x 10~ 5 



at 2 = 0.2467 and 
at 2 = 0.6847 



(3) 



assuming (perhaps incorrectly) a constant 
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4 DISCUSSION 

Very recently, COO has reported upper limits on Ay/y us- 
ing the same two absorption systems we analyse above. 
They obtained new H 1 21 cm spectra of each absorber and 
compared the measured redshift with the published molec- 
ular redshifts. COO find Ay/y = (+1.0 ± 0.3) x 10~ 5 for 
TXS 0218+357 and Ay/y = (+1.29 ± 0.08) x 10~ 5 for 
PKS 1413+135 where the errors are statistical only. They 
add to this an additional error of ±1 x 10~ 5 , which is de- 
rived from other measurement uncertainties in the H I data 
(e.g. frequency calibration uncertainties). However, COO ar- 
gue that the line of sight problem above could lead to errors 
as large as lOkms -1 on the basis of typical sub-kiloparsec 
ISM motions. They therefore conclude with an upper limit 
of \a/a\ < 3.5 x 10~ 15 yr~ 1 (assuming constant g p ) for a 
look-back time (within their assumed cosmology) of 4.8 Gyr 
to z = 0.6847. This corresponds to \Ay/y\ < 1.7 x 10" 5 . 
This final dominant error term of lOkms -1 is much greater 
than the 1.2 kms -1 error we apply in equations |^ and ^ but 
we stress that our value has been obtained empirically and 
so should be reliable. 

As mentioned in Section [l], tight constraints on Aa/a 
come from optical absorption line studies. One method is 
to compare the spacing of single alkali doublets (ADs) in 
QSO and laboratory spectra (Cowie & Songaila 1995; Var- 
shalovich, Panchuk & Ivanchik 1996). The most precise con- 
straint to date using this method was obtained by Murphy 
et al. (2001c), Aa/a = (-0.5 ± 1.3) x 10~ 5 , where we anal- 
ysed 21 Si IV doublets (2 < z < 3) in 13 QSO spectra (see 
Fig. 3). However, this method does not take advantage of 
the large relativistic shifts in the ground state energy since 
it compares transitions from the same ground state. 

Recently, a new many-multiplet (MM) method was sug- 
gested (Dzuba et al. 1999a,b; Webb et al. 1999) which al- 
lows an order of magnitude gain in precision over the AD 
method. The magnitude of the relativistic shifts is sensitive 
to the nuclear charge and so comparison of transitions in 
light and heavy ions yields tight constraints on Aa/a. Also, 
since any transition can be used, comparison of transitions 
with different ground states is possible. A further increase 
in precision comes from increased statistics since one is not 
restricted to using just a single AD. 

The first application of the MM method (Webb et al. 
1999) used 30 Mgl/Mgn/Fen absorption systems (0.5 < 
2 < 1.6) towards 17 QSOs and suggested that a may have 
been smaller in the past: Aa/a = (-1.09 ± 0.36) x 10~ 5 . 
More recently, we have re-analysed and confirmed the Webb 
et al. results with improved techniques and extended our 
analysis to a higher redshift sample of damped Lyman-a sys- 
tems (Murphy et al. 2001a; Webb et al. 2001). We compare 
our new mm/Hi results with these optical results in Fig. 3 
and we see that they are consistent. The weighted mean for 
the optical results, Aa/a = (-0.72 + 0.18) x 10" 5 , is signifi- 
cant at the 4.1a level and so a thorough search for systematic 
errors has been carried out (Murphy et al. 2001b). However, 
no effect was found that could explain the results. Fig. 3 
clearly shows that, assuming a constant g p , a large sample 
of mm/Hi absorption systems, over a somewhat larger red- 
shift range, will allow us to confirm or rule out the change 
in a implied by the optical data. 
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Figure 3. Comparison of our new mm/Hi constraints (crosses, 
from equation ^) and our recent optical constraints. The star rep- 
resents the weighted mean of Aa/a using 21 Si IV alkali doublets 
from Murphy et al. (2000c) while the dots represent a 7x 7 binning 
of the 49 absorption systems analysed with the many-multiplet 
method in Murphy et al. (2000a). 
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